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Abstract 

We have measured the conductivity of high-mobility (001) Si metal-oxide-semiconductor field 
effect transistors (MOSFETs) over wide ranges of electron densities n = (1.8 — 15) x 10 11 cm -2 , 
temperatures T = 30 mK - 4.2 K, and in-plane magnetic fields Bn = (0 — 5) T. The experimental 
data have been analyzed using the theory of interaction effects in the conductivity a of disordered 
2D systems. The parameters essential for comparison with the theory, such as the intervalley 
scattering time and valley splitting, have been measured or evaluated in independent experiments. 
The observed behavior of <r, including its quasi-linear increase with decreasing T down to ~ 0.4 K 
and its downturn at lower temperatures, is in agreement with the theory. The values of the Fermi- 
liquid parameter obtained from the comparison agree with the corresponding values extracted 
from the analysis of Shubnikov-de Haas oscillations based on the theory of magnetooscillations in 
interacting 2D systems. 

PACS numbers: PACS: 
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I. INTRODUCTION 



The beginning of the 80's witnessed a triumph of the one-parameter scaling theory of 
localization^ 2 - and the theory of electron-electron interactions in disordered conductors (for a 

n 

review, see Ref. |3|). As a result, the peculiar low-temperature behavior of the conductivity of 
numerous low-dimensional systems has been successfully attributed to quantum interference 
effects (see, e.g., Ref. I4). Curiously, application of these ideas to the two-dimensional (2D) 
electron liquid in Si MOSFETs - one of the most ubiquitous 2D systems - remained a chal- 
lenge for more than 25 years. In the early experiments^ with low-mobility (fj, ~ 0.1 m 2 /Vs) 
structures in the regime of relatively high electron densities (n > 10 12 cm~ 2 ), a decrease of 
the conductivity with cooling has been observed, apparently in qualitative agreement with 
the ideas of weak localization^ and electron-electron interactions^. However, the quantita- 
tive description of this behavior remained a problem^. The disagreement with the theory 
became qualitative with the advent of high-mobility (fi > 2 m 2 /Vs) structures: the low— T 
conductivity of high-/i Si MOSFETs increased with decreasing temperature^, in striking 
contrast to the behavior of many other 2D systems. This "metallic" behavior of the con- 
ductivity is especially pronounced at low electron densities n ~ (1 — 2) x 10 11 cm~ 2 , where a 
five-fold increase of cx(T) with cooling was observecr^^ 1 ^ 1 ^. Later, the quasi-linear increase 
of cr(T) was observed in many high- mobility 2D systems in the "dilute" regime, including 
GaAs heterostructures (both p-type^^^ and n-type^ 1 ^), Si/SiGe2£, and AlAs2i quan- 
tum wells. 

The subsequent development of theory and experiment led to significant progress in our 
understanding of the low-temperature transport in high-mobility systems in the regime of 
low electron densities. In the mid-80's, the quasi-linear "metallic" dependence observed in 
the ballistic interaction regime Tt ^> 1 (r is the transport mean free time, here and below we 
set h = ks = e = 1) was attributed^ to weakening of the screening of the scattering potential 
with increasing T. Observation of a strong negative magnetoconductance induced by the 
in-plane magnetic field 2 ^ and renormalization of the effective electron mass and g— factor in 
these structures" also suggested that the electron-electron interactions play an important 
role in this phenomenon. More recently, Zala, Narozhny, and Aleiner (ZNA}2£ developed a 
theory that took into account all interaction contributions to the conductivity, including the 
exchange ones. This theory offers a unified approach to both ballistic (Tt ^> 1) and diffusive 
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(TV <C 1) interaction regimes by considering the quantum interference between electron 
waves scattered off a short-range random potential "dressed" by Friedel oscillations of the 
electron density The theory was extended for the case of a long-range scattering potential 
by Gornyi and Mirlin (GM)2I. The theories^ 1 ^ naturally incorporate the Altshuler-Aronov 
results for the interaction corrections to the conductivity in the diffusive regime^. 

For the diffusive regime, a more general approach to interacting systems based on the 
nonlinear a-model has been developed by Finkel'stein^. Recently, the renormalization group 
(RG) equations of this theory^ 1 ^ (obtained in the first order in — and in all orders in 
interaction) have been compared with the conductivity of Si MOSFETs at low electron 
densities^^^^. 

The RG equations^'22'2&2^i describe the length scale (temperature) evolution of the 
resistivity and interaction parameters for a 2D electron system in the diffusive regime^. 
However, at high electron densities, the temperature range corresponding to the diffusive 
regime shrinks. In contrast, the theory of interaction corrections^ 1 ^ is applicable over a 
wider T range (that includes both ballistic and diffusive regimes) provided a 3> 1 and 
Aa/a <C 1; these assumptions are well justified at high densities. 

The theories^ 1 ^ predict that the magnitude and sign of the the interaction correction 
Aa(T, B) is determined by the value of the Fermi-liquid parameter Fq (which can be found 
by measuring the Shubnikov-de Haas (SdH) oscillations in weak magnetic fields perpendic- 
ular to the plane of a 2D structure^ 1 ^ or the magnetoresistance in strong in-plane magnetic 
fields^). In particular, it is expected that the cx(T) dependence becomes "metallic" when 
Fq is negative and its absolute value exceeds a certain threshold (see Sec. II). 

The experimental studies of the conductivity in various low-carrier-density 2D systems 
in the ballistic (high-temperature) regime are in agreement with the ZNA and GM theories. 
The "metallicity" in all these systems is enhanced at low n due to an increase of the absolute 



value of Fq (see, e.g., Refs. l37l]38l]39l ). In Si MOSFETs, the interaction effects are especially 
strong due to the presence of two nearly degenerate valleys in the electron spectrum^. 
This enhancement, however, diminishes if the temperature T becomes smaller than the 
valley splitting Ay and intervalley scattering rate t~ 1 . As a result, with lowering T, the 
"metallic" dependence of cr(T) is expected to become weaker or even to be replaced with 
an "insulating" one. To the best of our knowledge, this behavior has not been observed for 
Si MOSFETs prior to our work. 
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This paper aims to study in detail the conductivity of 2D electron liquid in high-mobility 
(001) Si MOSFETs over a wide temperature range (T = 0.03 — 4.2 K) that includes both the 
diffusive and ballistic regimes. In particular, we observed for the first time that the "metallic" 
increase of a with cooling is followed by the downturn of cr(T) at lower temperatures. For the 
purpose of comparison with the ZNA theory 2 ^, we studied the range of not-too-low densities, 
n = (1.8 — 15) x 10 11 cm 2 , where the temperature and magnetic field dependences Aa(T, B) 
can still be treated as small corrections to the Drude conductivity a p. In principle, no fitting 
parameters are required for comparison with the theory, because we have measured Fq, Ay, 
and t v in independent experiments. However, below we take a slightly different approach: 
we obtain the Fq(ti) values from fitting the Aa(T, B\\) dependences with the ZNA theory 2 ^, 
and show that these values are consistent with the corresponding values extracted from the 
analysis of SdH oscillations^. We have also revealed shortcomings of earlier analysis of 
Aa(T), reanalyzed the available data, and compared the extracted values of F^in) with 
corresponding values from other measurements. We conclude that the experimental data 
are well described by the theory of interaction corrections 2 ^ at intermediate temperatures 
T ~ (0.3 — 4.2) K. For a quantitative analysis at ultra-low temperatures (T < 0.3 K), the 
interaction correction theory should be modified by taking into account finite intervalley 
scattering rates. 

The paper is organized as follows. In Section II we briefly summarize the theoretical 
results 2 ^ for the interaction corrections to the conductivity of a two-valley system. The 
experimental data are presented in Section III, along with the data analysis and discussion. 
The summary is given in Section IV. 

II. INTERACTION CORRECTIONS TO THE CONDUCTIVITY 

A. Temperature dependence of the conductivity in zero magnetic field 

In the ZNA theory 2 ^, the corrections Aa ee to the Drude conductivity <jd = nr/nib (m& ~ 
0.205m e is the electron band mass in (001) Si MOSFETs; for more detail, see discussion in 
section IIII A 31 and references therein) were calculated in both ballistic and diffusive regimes 
for all orders of the interaction strength and the leading order in l/(EpT) and T/Ep. In 
particular, the theory reproduces the Altshuler-Aronov correction^ to the conductivity in the 
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diffusive regime. To adapt the theoretical results^ to the case of (001) Si MOSFETs, one 
should take into account that the electron spectrum in this system has two almost degenerate 
valleys^. In zero magnetic field, Aa ee (T) for a system with two degenerate valleys in the 
absence of intervalley scattering can be written as follows^: 

Aa ee (T) = 5a c (T) + l55a T (T). (1) 

Here 5o~c is the so-called "charge" contribution which combines Fock correction and the 
singlet part of Hartree correction, and 5<jt is the "triplet" contribution due to the triplet 
part of Hartree term. The valley index can be considered as a pseudo-spin in multi-valley 
systems^, and the valley degeneracy determines the number of triplet terms due to the spin 
exchange processes between electrons in different valleys. For a system with two degenerate 
valleys, the total number of interaction channels is 4 x 4 = 16, among them 1 singlet and 
15 triplet terms (for comparison, there are 1 singlet and 3 triplet terms for a single-valley 
system) . 

Below we assume that the scattering potential is short-ranged which is relevant to Si 



MOSFETs. According to Ref. |26|, the charge term does not depend on the details of inter- 
actions: 

whereas the magnitude and sign of the triplet term is controlled by the Fermi-liquid 
parameter— Ffi: 

The functions / and t in Eqs. (I2 f3]l describe the crossover between the diffusive ( Aa ee oc In T) 
and ballistic (Acr ee oc T ) regimes; outside the crossover region, they change the value of 
Aa Pf , by only a few percent. The explicit expressions for these functions can be found in 
Ref. |26j. The diffusive-ballistic crossover is expected over some temperature range near 



1 -I- F a 

T* = LLlo_. (4) 

2%r V ; 

Equations (T2"f3"j) describe the quantum corrections in a system with the conductance a ^> 
1/2% at temperatures well below the Fermi energy [T <C (1 + F$) Ep\ ■ 

The sign and the magnitude of Aa ee is controlled by the Fermi-liquid parameter Fq . For 
a rough estimate, deeply in the ballistic regime the InT terms and the crossover functions 



t and / in Eqs. (121131) can be omitted. For example, at Tt = 10, the functions / and t 
contribute ~ 4% and ~ 11%, respectively, to the linear-in-Tr ballistic terms of Eqs. fl2|3l) . 
By neglecting these terms, we find that for a single-valley system the linear dependence 
Ao- ee (T) in the ballistic regime becomes "metallic" (da/dT < 0) at 3F CT /(1 + F CT ) < -1 
or Fq < —0.25, whereas for a system with two degenerate valleys, the "metallic" Aa ee (T) 
dependences are expected for 15^/(1 + Fq) < —1 or Fq < —0.06. Thus, the valley 
degeneracy extends the range of Fq{ti), and hence, the range of carrier densities n where 
the conduction exhibits the "metallic" behavior. 

B. Magnetoconductivity in the in-plane magnetic field 

The in-plane magnetic field, being coupled mostly to electron spins, provides a useful tool 
for exploring the interaction effects in the low-temperature conductivity of Si MOSFETs^. 
When the Zeeman energy Ez = gb^BB (gb = 2 is the bare g-factor, /ig is the Bohr magneton) 
becomes much greater than T, the number of triplet terms that contribute to Aa ee (T) is 
reduced from 15 to 7. Similar reduction of triplet terms is expected for a valley splitting 
A v > T. These two effects have been accounted by the theory of interaction corrections^ 1 ^; 
in the presence of the magnetic field and/or valley splitting the interaction correction to the 
conductivity can be expressed as follows^: 

Aa ee (T,T,FZ,B h A v ) = Aa ee (T) + 2Aa z (E Z ,T) 

+ 2Aa z (A v , T) + Aa z {E z + A v , T) + Aa z (E z - A v , T), (5) 

where Aa ee (T) is given by Eq. ([T|). All the terms Aa z (Z, T) have a form 



Aa z (Z, T) = a(Z, T) - a(0, T) = 5a b (Z) + 5a d (Z) 
1 



7T 



kA Z 



2txT' 



f n 



+ m (Zt, Tt; F° 



(6) 



if the relevant energies Z Ep (Z stands for E z , Ay, and combinations E z ± Ay). The 
explicit expressions for the functions Kb and Kd are given in Ref. [26j. In particular, Eq. ([6]) 
describes the interaction-driven magnetoconductivity in the magnetic fields which are much 
weaker than the field of full spin polarization of a system. Below we will neglect the function 
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m(0, Tr; Fq) which describes the crossover between the ballistic and diffusive regimes: this 
function is numerically small and does not modify the value of Aa(Z, T) outside the ballistic- 
diffusive crossover region by more then one percent. 

It is worth mentioning that in the framework of the RG theory, the magnetoconductance 
can also be described by the Castellani-Di Castro-Lee formular^i which is equivalent to 
Eq. (J5]) in the diffusive limit at Ay = 0. However, for the analysis of our magnetoconductivity 
data measured over a wide T-range that includes both diffusive and ballistic regimes, the 
interaction correction theory^ is more appropriate than the RG theory^. 

The interaction correction theory^ (as well as the RG theory^ 1 ^) does not take into 
account intervalley scattering. This approximation is valid when the intervalley scattering 
rate t~ 1 is much smaller than T. In the low-temperature limit T -C r" 1 , the electron states 
in different valleys are completely intermixed at the time scale ~ T _1 and the correction 
Aa ee (T) for a two- valley system is expected to coincide with that for a single- valley system. 

Since the interaction corrections to the conductivity Aa ee (T, B) depend on several param- 
eters such as r, Fq, Ay and t v , for testing the theoretical results, it is crucial to determine 
these parameters in independent measurements. This program is realized in the following 
section. 

III. SAMPLE CHARACTERIZATION AND DATA ANALYSIS 
A. Sample Characterization 

We have studied the temperature and magnetic field dependences of the conductivity 
for high-mobility (001) Si MOSFETs, which demonstrate the "metallic" quasi-linear cr(T) 
dependences at intermediate temperatures over a wide range of electron densities n. In this 
paper we present the data for two (001) Si MOSFETs, Si6-14 and Sil-46, with the gate oxide 
thickness 190±20nm and peak mobility //(0.1K) ~ 2m 2 /Vs. The a(T, B) dependences were 
measured over the temperature range T = 0.03 — 4.2 K using a standard low-frequency four- 
terminal technique. The measuring current was chosen sufficiently small (J ~ 1 — 3nA) 
to avoid overheating of electrons within this temperature ranged Our experimental set- 
up allowed us to independently control the magnetic field normal to the plane of 2D layer 
(B± ~ — 1.5. .. + 1.5 T) and the in-plane magnetic field (B\\ ~ — 8... + 8T); this cross-magnetic 
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field technique has been described in Ref. |43|. Unless otherwise stated, an in-plane magnetic 
field of Bn = 0.02 T was applied to quench the superconductivity in the current /voltage 



contact pads and the gate electrode which are made of thin aluminum films. 



1. Effective mass m* and g— factor. Interaction corrections to the magnetooscillations. 



The electron density n, the electron temperature T, the effective electron mass m*, and 
g— factor g* have been found from measurements of SdH oscillations (see also Ref.l25l). For 
fitting the oscillations, as the first step of the analysis, we have used Lifshitz-Kosevich (LK) 
formula^ which is valid for non-interacting 2D electrons if the amplitude of oscillations is 

small^: 

8a xx (u c ,T) ^ LK 
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(7) 



uj c J sinh(2ir 2 sT/uj c ) 

Here u c = Bj_/m* is the cyclotron frequency, Tjj = 1/2tit£) is the Dingle temperature, td 
is the elastic quantum scattering time. Figure 1 shows that the temperature dependence 
of the amplitude of the first harmonic, A\(T), is in agreement with Eq. (J7J) down to the 
lowest temperatures; this indicates that the electrons are not overheated (with respect to 
the thermal bath) by the bias current and/or noise. 

The Fq(t s ) values can be found from the measurements of the renormalized g-factor: 
Fq = (gb/g*) — 1- Here r s = l/W^y/wri is the dimensionless ratio of the Wigner-Seitz radius 
to the effective Bohr radius, a* B = ~K/me 2 , k = 7.7 is the average dielectric constant of Si 
and SiC>2, and m = 0.1 9m e is the electron band mass in bulk Si^ 1 ^. In our experiments^, 
g*(n) was obtained from the analysis of SdH oscillations as the ratio of two quantities: the 
measured renormalized electron spin susceptibility x*/Xb = g*m* ' / gbf^b and the effective 
mass m* (xb is the bare value of spin susceptibility). Observation of the beatings of SdH 
oscillations in crossed magnetic fields offers a straightforward (model-independent) method 
of finding g*m* On the other hand, an estimate of m* is based on a model-dependent 
analysis of the damping factor for the first harmonic of SdH oscillations, Ai(T, = const). 
According to the LK theory, the damping factor can be expressed as 



In [A\ K {T,B A 
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(8) 
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Our experiments show that In A\(T,B±) varies linearly with T within the experimental range 



T = (0.4 — 0.8) K (see Fig. 3 in Ref.l25l); this however does not prove the applicability of the 
LK theory, which disregards the interaction effects. It was recently shown^^Z, that due to 
the interference between electron-electron and electron-impurity interactions, the damping 
factor acquires an additional term in both the diffusive and ballistic regimes 

- In [Ai(T, B ± )] = (T + T D ) - a(T), (9) 



where 



and 



27r 2 m* 



tt (T) = -T^-T D (^-^), (10) 
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A=(l + ^-)^— (11) 

for a system with two degenerate valleys. 

The equation for Sm*(T)/m* resembles the one-loop renormalization of the effective mass 
(or Z) in the RG theory^ 1 ^ 1 ^. Our numerical simulations show that within the relevant 
interval T = (0.03 — 0.8) K and r s < 6, the In T terms in Eqs. (iTTj) can be replaced with a T- 
independent constant. By combining the LK result with the interaction- induced corrections 
and replacing all terms oc InT by a constant within our limited T range, we obtain the 
following linearized equation in the ballistic regime for the short-range scattering (td ~ t): 

- In [Ax(T, B ± )] =T + T D (1- 2nATr) =T + T D (l- i^^) . (12) 

In this case, the T— dependent correction to the Dingle temperature, <5T j^ T - > , is one half 
of the interaction correction to the conductivity^ (this factor | originates from the 

difference between the interaction corrections to the momentum relaxation and quantum 
scattering times^) . We note that the empirical procedure used for finding m* in our earlier 



paper (Ref.l25l) was based on the assumption that T£, = Tp [1 — 5a(T)/aD} 1 which differs 
from Eq. (112j) by a factor of \. 

At relatively high densities (which correspond to r s < 4), the corrections to the LK result 
are insignificant within the studied T-range. As r s increases, the temperature dependences 



of the oscillation magnitude predicted by the LK theory Eq. (jSj) and the interaction theory^ 
start deviating from each other. The values of \Fq\ extracted from SdH data using Eq. ffT2]) 
are larger than those obtained with the LK theory but smaller than \Fq\ obtained with 



the empirical procedure of Ref. 



25 



(e.g., at r 



6.2, the values of Fq obtained according to 



Eq. [Inland the empirical procedure of Ref. |25| are —0.40 and —0.45, respectively). We have 



reanalyzed the data of Ref.[25| using Eq. (Tl2|) and compared the corresponding values with 



Fq{ji) extracted from the a(T,B\\) dependences using the ZNA theory (see below). 
2. Valley splitting and intervalley scattering 

The analysis of SdH oscillations using Eq. ([7]) also allowed us to estimate the energy 
splitting Ay between the valleys. A non-zero valley splitting results in the beating of SdH 
oscillations^. Figure 1 (a,b) shows the SdH oscillations for samples Si6-14 and Sil-46 (the 
electron densities are 6.1 x 10 n cm~ 2 and 1 x 10 12 cm~ 2 , respectively). The amplitude of 
SdH oscillations normalized by the first harmonic A\ is expected to be field-independent if 
Ay = 0. A noticeable reduction of the SdH amplitude observed for both samples at small 
fields can be attributed to a finite valley splitting. Although the node of SdH oscillations 
(expected at B « 0.15 T) cannot be resolved for samples with mobilities fi = 2m 2 /Vs, 
Ay can be estimated from fitting of the B— dependence of the SdH amplitude with Eq. ([7]) 
modified for the case of a finite Ay: Ay = 0.4 K for sample Si6-14 and 0.7 K for Sil-46. 
This estimate provides the upper limit for Ay at B — 0: in nonzero B± fields, Ay may be 
enhanced by the inter-level interaction effects^ 1 ^. 



The intervalley scattering rate for sample Si6-14 was measured earlier in Ref.l52l by ana- 
lyzing the weak-localization (WL) magnetoresistance. It was found that t v is temperature- 
independent and the ratio t v /t decreases monotonically with increasing electron density. 
For Si6-14 at n = (3 — 6) x 10 11 cm -2 , t v = 20 ps (r" 1 = 0.36 K) is approximately ten times 
greater than the transport time r = 2 ps. 

3. Relaxation time r and the band mass mb 

The momentum relaxation time r was determined from the Drude conductivity gd = 
nr/mi,, which was found by extrapolating the quasi-linear a(T) dependence observed in 
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the ballistic regime to T = 0. Note that in order to extract r from the Drude conductiv- 
ity, one should use the bare mass m&: according to the Kohn theorem, the response of a 
translationally-invariant system to the electromagnetic field is described by mj in the pres- 
ence of electron-electron interactions; this result also holds for weak disorder (EpT ^> 1). It 
is worth mentioning that several prior publications^ 1 ^ 1 ^, including our paper—, incorrectly 
used m* instead of mj to estimate r from gd\ this affects the value of the fitting parameters 
extracted from comparison with the theory 2 ^ as shown below. 

The textbook valued for the light electron mass in bulk Si is m!f D ^ « 0.19m e . For 
inversion layers on (001) Si-surface, Kunze and Lautz^ have obtained m[ 2D ^/m e = (0.19 — 
0.22)±0.02 from tunneling measurements. Our recent m*(n) data obtained from the analysis 
of SdH oscillations^ over a wide range of densities r s = 1.4 — 8.5, can be fitted with a 
polynomial m*(r s ) = 0.205m e (l + 0.035r s + 0.00025rg). These m*/m e data agree well with 
earlier values of m* extracted from SdH oscillations^^ 1 ^ 1 ^ in narrower ranges of densities. 
By extrapolating the polynomial m*(r s ) to r s = we obtain mf D) /m e = 0.205 ± 0.005, the 
value which we adopted throughout the paper™; available measurements of the cyclotron 
resonance^ do not contradict and do not refine this value. 

In principle, the aforementioned complete characterization of samples allows us to com- 
pare the Aa(T, B) dependences with the theory 2 ^ without any fitting parameters (with a 
caveat that the theory does not take into account the intervalley scattering, see the discus- 
sion below). However, throughout this paper we adopt an equivalent, but more convenient 
procedure: for each electron density, Fq{ji) will be considered as a single parameter for 
fitting Aa(T,B), and these values of Fq will be compared with the corresponding values 
obtained from the SdH oscillations^. 

B. Temperature dependences of the conductivity at Bu = 

The temperature dependences of the conductivity cr(T) for sample Si6-14 are shown in 
Fig. 2. In these measurements, we applied a fixed B± = 0.1 T that is sufficient to suppress the 
temperature dependence of the WL correction in the studied temperature range. The cr(T) 
dependences are non-monotonic for all studied densities [n = (1.8 — 15) x 10 11 cm -2 for Si6-14 
and n = (10 — 15) x 10 n cm -2 for Sil-46]: a quasi-linear increase of a with cooling, observed 
down to ~ 0.5 K, is replaced at lower T with a decrease of a. Note that in our previous 
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experiments^, we observed a trend of cr(T) saturation at T < 0.4 K rather than the decrease 
of the conductivity. One of the reasons for this might have been "heating" of electrons by 
high frequency noise: only after thorough filtering of all leads connected to the sample were 
we able to decrease the electron temperature down to ~ 30mK. Similar downturn of a(T), 
although at much lower temperatures, has been recently observed in high-/i GaAs FETs 
at low electron densities^. We note also that earlier, a downturn of cr(T) was observed in 
Si-MOSFETs^i but at much higher electron densities (n > 30 x 10 11 cm -2 , r s < 1.4) and 
at much higher temperatures T ~ 10K (Tr rs 1). For such high densities the interaction 
corrections to a become negative and the downturn of a(T) was related to the crossover 
between a "metallic" high-temperature a(T) dependence (that is due to electron-phonon and 
intersubband scattering effects) and resulting "insulating" low-temperature a (T) dependence 
(that is due to weak localization and negative interaction corrections contributions to a(T)). 

Below we use the following strategy for analyzing the Aa(T) dependences. First, we find 
Fq by fitting the quasi-linear cr(T) dependences observed in the ballistic regime (T > 0.5 K) 
with Eqs. dUE]). The effect of valley splitting and intervalley scattering on Aa ee (T) can be 
neglected at T > A^r" 1 and the analysis is significantly simplified. The corresponding 
values of Fq(ti) are shown in Fig. 5. The T range available for fitting in this regime "shrinks" 
rapidly at low n: the growth of \Fq\ and decrease of Ep with decreasing n, lead to violation of 
the condition T < (l+F°) 2 E F (e.g., at n = 1.8 x 10 cm this occurs at temperatures above 
2K). This might be one of the reasons for the observed deviation of the high-temperature 
a(T) from the linear-in- (T/Ep) theory^. Also, the higher-order corrections might become 
significant at low n when Aa(T)/ao ~ 1 (see Fig. 2). 

After finding the Fq values (which are temperature-independent in the studied tempera- 
ture range), we proceed with the analysis of the low— T part of the cr(T) dependences, where 
the crossover from da/dT < to da/dT > was observed. We note that the crossover occurs 
when the temperature becomes smaller than two characteristic temperature scales - Ay and 
r" 1 - which are of the same order of magnitude for the studied structures. We emphasize 
that according to the ZNA theory^, the ballistic-diffusive crossover should not lead to the 
change of the sign of da/dT. In contrast, the valley splitting and the intervalley scattering 
may result in the sign change for da/dT because these processes reduce the number of triplet 
components at T < Ay and T < r" 1 . 

The theory^ takes into account a finite Ay but not r^ 1 . The solid red curves in Fig. 2 
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are calculated for Ay = 0.4 K (the estimated value of Ay for sample Si6-14) and r" 1 = 0. 
It is clear that the change in the number of triplet components from 15 (T ^> Ay) to 7 
(T <C A„)2£i£2 ( S ee also Eq. ED is not sufficient to explain the shape of the cr(T) downturn. 
The effect of strong intervalley scattering is illustrated in Fig. 2 by dashed green curves 
calculated with 3 triplet components (to model roughly the case of T <C t" 1 when the 
valleys are completely intermixed). In the absence of a detailed theory that would account 
for intervalley mixing, we attempted to fit the experimental data with an empirical crossover 
function for the number of triplet components N trip \ ct (x) = 9 + 6[exp(— 0.3/x) — exp(— 30x)], 
where x = Tt v . This crossover function provides correct asymptotic limits for Applet : 3 at 
T <C r" 1 and 15 at T ^> r" 1 . Figure 2 shows that using this function, we can reasonably 
well describe the shape of experimental cx(T) dependences for all studied electron densities. 

C. Temperature dependences of the conductivity at non-zero B» 

Better understanding of different contributions to Acx(T) can be achieved by measuring 
the conductivity in strong in-plane magnetic fields B\\ ^> T j ' g^B- The evolution of exper- 
imental dependences c(T) with Bu is shown in Fig. 3 for two samples at different electron 
densities. The theoretical curves in Fig. 3 were calculated using the Fq(ti) values extracted 
from the analysis of 5a ee (T, B\\ = 0) (see Fig. 2). The transport time r was calculated for 
each B\\ value from the Drude conductivity <td (J3n) which in turn was estimated by extrap- 
olating the quasi-linear part of the a(T, B») dependence to T = 0. The observed behavior is 
in line with our analysis of the cr(T) dependences in Section III B. Indeed, the magnitude of 
the triplet contribution is expected to be reduced when the Zeeman energy becomes greater 
than T. This effect is more pronounced within the range Ay, t" 1 < T < g^fisB, where a 
strong magnetic field reduces the number of triplet components from 15 to 7. For example, 
at n = 1 x 10 12 cm -2 (see Fig. 3 d) the "metallic" behavior disappears at T < 1 K and 
B\\ = 3 T, which is in agreement with the theory. At lower T, the number of triplet com- 
ponents is smaller than 15 even at B» = due to valley splitting and intervalley scattering, 
and the effect of B\\ on Aa(T) is less prominent. 
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D. Magnetoconductivity 



To test the theoretical predictions on the magnetoconductivity (MC) induced by in-plane 
magnetic fields, we also measured the cr(B\\) dependences at fixed T. Similar measurements 



have been performed in the past (see, e.g., Refs. l37lJ38li63l ). but no detailed comparison with 
the theory was carried out at that time. The MC for sample Si6-14 over the field range 
—4.5 < JBm < 4.5 T is shown for different densities and temperatures in Fig. 4. In these mea- 
surements, special care was taken to reduce the magnetic field component perpendicular to 
the plane of the structure: even a 1° misalignment between the sample's plane and the mag- 
net axis (which results in B± ~ 50 G at B» = 3 T) may be sufficiently strong for suppressing 
the WL corrections at low T. To eliminate B±, we used the cross- magnetic-field set-up^. 
For each value of B\\, we measured the dependence cr(B±) by sweeping B± and recorded 
the minimum value of cr(B±) which corresponded to the zero WL magnetoconductance and, 
thus, B± = 0. This method allowed us to compensate B± with an accuracy better than 
10 G. 

The theoretical Aa(B») dependences (see Eqs. (15fB]) ). plotted in Fig. 4 as solid curves, 
describe the observed MC very well in not-too-strong magnetic fields g^BB\\ < 0.2E F . 
Again, as in the case of fitting the Acr(T) dependences, the only adjustable parameter 
was the Fq{ji) value extracted for each density from fitting the MC at high temperatures 
(« 0.7 K) where the effects of valley splitting or intervalley scattering on Aa ee (T, B) can 
be neglected. Note that all the theoretical curves plotted in Fig. 4 for the same n were 
calculated for a fixed Fq(ti), i.e. neglecting possible dependence Fq(B»). The detailed 



analysis of the spin susceptibility x* ^ 9*™* in strong magnetic field, presented in Ref. 
shows that the product g*m* decreases with an increase of B\\ by as much as ~ 20%. Our 
estimate shows that by ignoring the g*(B) dependence, we might reduce the value of \Fq\ 
by ~ 10% (see below), which is close to the accuracy of extraction of F£ from the data 
in strong magnetic fields. As B\\ grows and/or n decreases, the data start deviating from 
the theoretical curves (see Fig. 4d); this deviation can be attributed to the violation of the 
condition g^sB <C E f required for applicability of Eqs. (15f5]) . 
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E. The Fq(h) dependence 

The Fq values obtained from fitting the Aer(T) and Aa(B») dependences with the 
theory^ are shown in Fig. 5. For comparison, we have also plotted the Fq values obtained 
from the analysis of SdH oscillations measured for sample Si6-14 using the theories^ 1 ^ (see 
Section III A). The Fq values obtained from fitting Acx(T) are in good agreement with the 
corresponding values extracted from the analysis of SdH data. 

At the same time, the \Ftf\ values obtained from fitting the Aa(B\\) dependences at r s > 4 
are systematically smaller than the corresponding values obtained from fitting Acx(T) and 



37 



38 



There are at least two factors 



SdH oscillations. This trend was earlier reported in Refs. 
that can reduce this discrepancy. One of them, a potential decrease of g* in strong B\\, was 
mentioned in Section III D. The other factor is more subtle. In our analysis, we neglected 
the dependence of the WL correction <5<twl on B\\. However, our measurements show that 
5<twl decreases with an increase of the in-plane magnetic field, which leads to a positive 
magnetoconductance. There are at least two potential reasons for this dependence: (a) the 
Si-Si02 interface roughness transforms a uniform in-plane field into a random perpendicular 



field (see, e.g., Refs. 65,66 IJ67J and references therein), and (b) a finite extent of electron wave 



functions in the direction perpendicular to the plane of a quantum well causes sub-band 



mixing by the magnetic field and disorder (see, e.g., Ref.|68| and references therein). 

Phenomenologically, both effects can be described in terms of a decrease of the dephasing 
length L v with B\\. For example, from the analysis of the WL magnetoresistance mea- 
sured for different values of B» for sample Si6-14 at n — 1 X 10 12 cm~ 2 and T = 0.3K, we 
have extracted L v (Bh = 0) = 1.3fim and L V (B» = 3T) = 0.8fim. Our estimates show 
that the positive magnetoconductance associated with the suppression of 5<twl by B\\ can 
account for ~ 10% of the observed .By-induced magnetoconductance. After taking the de- 
pendences g*(B\\) and 8a^fj J (Bn) into account, the \Fq\ values extracted from the MC should 
be increased by ~ 20%. The corresponding downshift of the Fq values extracted from the 
Aa(B) data is shown by arrows in Fig. 5 a, it significantly reduces the discrepancy between 
the values of Fq extracted from the B\\ = data and the data measured at B» — 1 — 5 T. 

Finally, in order to compare our data with other available Aa(T) data for high-mobility 
Si MOSFETs, we have used Eqs. (JTH6]) to estimate the Fq values from the quasi-linear 
<t(T, B = 0) dependences measured by Shashkin et al.— and Vitkalov et al.— . In this 
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analysis, we estimated r from a(T — > 0) using the band mass rather than m*. As for the 



va. 



38 



ley splitting, we have used Ay = 0.6 K and 1 K for the a na 
respectively (The value Ay ~ (1 — 2) K used in Ref. 



36 



ysis of data from Ref. 



36 



and 



seems to be too large, as it 



wou 
Ref. 



d lead to the appearance of beating of SdH oscillations in the field range studied in 



38 



691 ). We have also taken into account the WL correction neglected in both Ref. 



36 and 



. Figure 5 b shows that the Fq(t s ) values estimated for different Si MOSFETs using the 
ZNA theory are in good agreement with each other. 

The Fq values obtained on the basis of the interaction correction theory and plotted 
in Fig. 5 may be compared also with the values of Fq = —72/(1 + 72) predicted by the 
RG theory and measured experimentally in Refs. l3l|]32l . Extrapolation of our Fq(t s ) data 
to lower densities, provides the value Fq « —0.5 at n m 1 x 10 11 cm -2 (r s ~ 8), which 
is smaller than the value Fq = —0.31 (72 = 0.45) predicted by the one-loop RG theory 
for the temperature T max corresponding to the p(T) maximum^ 1 ^ 1 ^ (e.g., T max « 3K for 
n = 1.2 x 10 n cm- 2 ). 

The experimental test£^2 of the RG theory was conducted at temperatures higher than 
that in the experiments described in the present work. Within the framework of the RG 
theory, the interaction parameter 72 is expected to increases with decreasing j££i2ii22, i n 
principle, it can reach at T < IK the value of ~ 1 which corresponds to F$ = —0.5. (Note 
that the factor-of-two increase of 72 (from 0.45 to 1) is beyond the range of the applicability 
of the one- loop RG theory.) Another problem is that the spin susceptibility x* ^ 9*™* 
obtained from the SdH data is almost T-independentr^, in contrast to the expected increase 
of 72 (and, hence, \Fq\ and g*) with cooling. This contradiction can be resolved if the 
T-dependence of g* is exactly compensated by the opposite T-dependence of m*, so that 
X* oc g*m* remains almost constant. The reason for this compensation is not clear and 
requires both experimental and theoretical studies. 



IV. CONCLUSION 



Our experiments show that the low-T behavior of the conductivity of high-mobility 
(001) Si MOSFETs is well described by the theory of interaction effects in systems with 
short-range disorder—. Over a wide range of intermediate temperatures (Ay, r" 1 , gb^sB < 
T Ep), the interaction effects are strongly enhanced in Si MOSFETs due to the presence 
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of two valleys in the electron spectrum. This factor, in combination with the interaction- 
driven renormalization of the Fermi-liquid parameter F£, leads to an increase of a with 
decreasing T. At lower temperatures (T < A^r" 1 , g b fi B B, E F ), the triplet contribution 
to Acr ee (T) is significantly reduced due to valley splitting and/or intervalley scattering. As 
a result, the "metallic" behavior of a is replaced with a more conventional, "insulating" 
behavior. The Ftf values obtained from fitting the experimental data with the theory 2 ^ 
agree well with the Fq data obtained from the analysis of SdH oscillations in these samples. 
However, it remains unclear how to reconcile the Fq values obtained at low n from fitting 
the cr(T) and SdH data by using the interaction correction theory with the corresponding 
values obtained within framework of the RG theory. 

We emphasize that for the detailed analysis of the interaction-induced contributions to 
the conductivity, it is important to measure such parameters as the valley splitting and 
intervalley scattering rate in independent experiments. Finally, for a quantitative description 
of the interaction effects to the conductivity Aa ee at low temperatures, both the interaction 
correction theory and RG theory should be extended to the case of a finite intervalley 
scattering rate. 

V. ACKNOWLEDGMENTS 

The authors are thankful to E. Abrahams, I. Aleiner, I. Burmistrov, I. Gornyi, G. Kotliar, 
and A. Mirlin for illuminating discussions. The research at Rutgers University was partially 
supported by the NSF grant ECE-0608842. The research at Lebedev Institute was supported 
by RFBR, Programs of the RAS, Russian Ministry for Education and Science, and the 
Program "Leading Scientific Schools". 



1 E. Abrahams, P. W. Anderson, D. C. Licceardello, and T.V. Ramakrishnan, Phys. Rev. Lett. 
42, 673 (1979). 

2 F. J. Wegner, Phys. Rev. B 19, 783 (1979). 

3 B.L. Altshuler, A. G. Aronov, and PA. Lee, Phys. Rev. Lett. 44, 1288 (1980); B.L. Altshuler 
and A. G. Aronov, in Electron-Electron Interaction in Disordered Systems, edited by A. L. Efros 
and M. Pollak (North-Holland, Amsterdam, 1985). 

17 



4 B. L. Altshuler, A. G. Aronov, M. E. Gershenson, and Yu. V. Sharvin, Sov. Sci. Rev. A. Phys. 
9, 223-254 (1987). 

5 D. J. Bishop, D. C. Tsui, and R. C. Dynes, Phys. Rev. Lett. 44, 1153 (1980); Phys. Rev. B 26, 
773 (1982). 

6 V.T. Dolgopolov, S.I. Dorozhkin, A. A. Shashkin, Sol. State Commun. 50, 273 (1984). 

7 L.P. Gorkov, A.I. Larkin, and D.E. Khmelnitskii, JETP Lett., 30, 228 (1979). [Pis'ma v Zh. 
Eksp. Teor. Fiz. 30, 248 (1979)]. 

8 M.S. Burdis and C.C. Dean, Phys. Rev. B 38, 3269 (1988). 

9 K.M. Cham and R. G. Wheeler, Phys. Rev. Lett. 44, 1472 (1980). 

R.P. Smith and P.J. Stiles, Sol. State Commun. 58, 511 (1986). 

1 T. N. Zavaritskaya and E. I. Zavaritskaya, JETP Lett. 45, 609 (1987). [Pis'ma v Zh. Eksp. Teor. 
Fiz. 45, 476 (1987)]. 

2 E. A. Virodov, V.T. Dolgopolov, S.I. Dorozhkin, N.B. Zhitenev, Zh. Eksp. Teor. Fiz. 94, 234; 
[JETP 67, 998 (1988)]. 

3 S. V. Kravchenko, G. V. Kravchenko, J. E. Furneaux, V. M. Pudalov, and M. D'lorio, Phys. Rev. 
B 50, 8039 (1994). 

4 Y. Hanein, U. Meirav, D. Shahar, C. C. Li, D. C. Tsui, and H. Shtrikman, Phys. Rev. Lett. 80, 
1288 (1998); Y. Hanein, D. Shahar, J. Yoon, C. C. Li, D. C. Tsui, and H. Shtriktman, Phys. Rev. 
B 56, R7520 (1998); A. R. Hamilton, M.Y. Simmons, M. Pepper, E. H. Linfield, P.D. Rose, 
and D.A. Ritchie, Phys. Rev. Lett. 82, 1542 (1999); A. P. Mills, A. P. Ramirez, L.N. Pfeiffer, 
and K. West, Phys. Rev. Lett. 83, 2805 (1999). 

5 M.Y. Simmons, A. R. Hamilton, M. Pepper, E. H. Linfield, P.D. Rose, D.A. Ritchie, 
A.K. Savchenko, and T. G. Griffiths, Phys. Rev. Lett. 80, 1292, (1998); J. Yoon, C.C. Li, 
D. Shahar, D. C. Tsui, and M. Shayegan, Phys. Rev. Lett. 82, 1744 (1999). 

6 S.J. Papadakis, E.P. De Poortere, H.C. Manoharan, M. Shayegan, and R. Winkler, Science 
283, 2056 (1999). 

7 Y. Hanein, D. Shahar, H. Shtrikman, J. Yoon, C. C. Li, and D. C. Tsui, Nature 400, 735 (1999). 

8 Y. Hanein, D. Shahar, J. Yoon, C. C. Li, D. C. Tsui, and H. Shtrikman, Phys. Rev. B 58, 
R13338 (1998). 

9 E. Ribeiro, R. D. Jaggi, T. Heinzel, K. Ensslin, G. Medeiros-Ribeiro, and P. M. Petroff, Phys. 
Rev. Lett. 82, 996 (1999). 

18 



20 P.T. Coleridge, R.L. Williams, Y. Feng, and P. Zawadzki, Phys. Rev. B 56, R12764 (1997); 
J. Lam, M. D'lorio, D. Brown, and H. Lafontaine, Phys. Rev. B 56, R12741 (1997); P.T. Co- 
leridge, A. Sachrajda, and P. Zawadzki, Phys. Rev. B 65, 125328 (2002). 

21 S.J. Papadakis and M. Shayegan, Phys. Rev. B 57, R15068 (1998). 

22 F. Stern and W. E. Howard, Phys. Rev. 163, 816 (1963); A. Gold and V.T. Dolgopolov, Phys. 
Rev. B 33, 1076 (1986); S. Das Sarma, Phys. Rev. B 33, 5401 (1986); S. Das Sarma, and 
E.H. Hwang, Phys. Rev. Lett. 83, 164 (1999); S. Das Sarma, and E.H. Hwang, Phys. Rev. B 
69, 195305 (2004). 

23 D. Simonian S.V. Kravchenko, M. P. Sarachik, and V.M. Pudalov, Phys. Rev. Lett. 79, 2304 
(1997); V.M. Pudalov, G. Brunthaler, A. Prinz, G. Bauer, JETP Lett. 65, 933 (1997). 

24 T. Okamoto, K. Hosoya, S. Kawaji, and A. Yagi, Phys. Rev. Lett. 82, 3875 (1999). 

25 V. M. Pudalov, M. E. Gershenson, H. Kojima, N. Butch, E. M. Dizhur, G. Brunthaler, A. Prinz, 
and G. Bauer, Phys. Rev. Lett. 88, 196404 (2002). 

26 G. Zala, B.N. Narozhny, and I.L. Aleiner, Phys. Rev. B 64, 214204 (2001); ibid 65, 020201 
(2002). 

27 I.V. Gornyi and A.D. Mirlin, Phys. Rev. B. 69, 045313 (2004). 

28 A.M. Finkel'stein, Z. Phys. B 56 189 (1984); Sov. Sci. Rev. A, Phys. 14, 1-100, (1990). 

29 A. Punnoose and A.M. Finkel'stein, Phys. Rev. Lett. 88 16802 (2002). 

30 A. Punnoose and A.M. Finkel'stein, Science 310, 289 (2005). 

31 D.A. Knyazev, O.E. Omel'yanovskii, V.M. Pudalov, and IS. Burmistrov, JETP Lett. 84, 662 
(2006). 

32 S. Anissimova, S.V. Kravchenko, A. Punnoose, A.M. Finkel'stein, and T.M. Klapwijk, Nature 
Phys. 3, 707 (2007). 

33 D.A. Knyazev, O.E. Omel'yanovskii, V.M. Pudalov, I. S. Burmistrov, Phys. Rev. Lett. 100, 
046405 (2008). 

34 C. Castellani, C.Di Castro, P. A. Lee, and M. Ma, Phys. Rev. B 30, 527 (1984). 

35 C. Castellani and C.Di Castro, and P. A. Lee, Phys. Rev. B 57, 9381 (1998). 

36 A. A. Shashkin, S.V. Kravchenko, V.T. Dolgopolov, and T.M. Klapwijk, Phys. Rev. B 66, 
073303 (2002). 

37 V. M. Pudalov, M. E. Gershenson, H. Kojima, G. Brunthaler, A. Prinz, and G. Bauer, Phys. 
Rev. Lett. 91, 126403 (2003). 

19 



38 S.A. Vitkalov, K. James, B.N. Narozhny, M. P. Sarachik, and T.M. Klapwijk, Phys. Rev. B 
67, 113310 (2003). 

39 W.R. Clarke, C.E. Yasin, A. R. Hamilton, A. P. Micolich, M.Y. Simmons, K. Muraki, Y. Hi- 
rayama, M. Pepper and D. A. Ritchie, Nature Physics 4, 55 (2007). 

40 T. Ando, A.B. Fowler, and F. Stern, Rev. Mod. Phys. 54, 437 (1982). 

41 I. S. Burmistrov and N. M. Chtchelkatchev, JETP Lett. 84, 558 (2006); Phys. Rev. B 77, 195319 
(2008). 

42 O. Prus, M. Reznikov, U. Sivan, and V. M. Pudalov, Phys. Rev. Lett. 88, 016801 (2002). 

43 M. E. Gershenson, V. M. Pudalov, H. Kojima, N. Butch, G. Bauer, G. Brunthaler, and A. Prinz, 
Physica E 12, 585 (2002). 

44 I.M. Lifshitz, A.M. Kosevich, Zh. Eksp. Teor. Fiz. 29, 730 (1955); A. Isihara and L. Smrcka, 
J. Phys. C 19, 6777 (1986). 

45 G.W. Martin, D.L. Maslov, and M.Yu. Reizer, Phys. Rev. B 68, 241309(R) (2003). 

46 Here we used a traditional value of nib m order to fasilitate the comparison with the results of 
other experimental works. 

47 Y. Adamov, I.V. Gornyi, and A.D. Mirlin, Phys. Rev. B 73, 045426 (2006). 

8 Interaction correction to the conductivity [Eqs. (Q3I3D] in the ballistic regime has the form 
Sa(T)/a D = Tt[1 + 15F<f/(l + Fg)]/ircr D . 

49 The quantum correction to the transport scattering rate r" 1 differs from that for the quantum 
scattering rate r^ 1 by the term (1 — cos(cf))) in the integrand, where 4> is the scattering angle. 
According to ZNA theory 2 ^, the interaction corrections to the conductivity are determined by 
backscattering events for which (f> ~ ir, or (1 — cos{<p)) ~ 2. 

50 V. M. Pudalov, A. Punnoose, G. Brunthaler, A. Prinz, G. Bauer, |cond-mat / 0104347| In this 



paper, three Si MOSFET samples with anomalously large values Ay ~ 1 — 2K were studied; 
the high Ay values enabled observation of the SdH oscillation beats in low magnetic fileds. 

51 V. M. Pudalov, S. G. Semenchinskii, V. S. Edelman, Pis'ma v ZhETF, 41(6), 265 (1985). [JETP 
Lett. 41, 325 (1985)]. 

52 A. Yu. Kuntsevich, N.N. Klimov, S. A. Tarasenko, N. S. Averkiev, V. M. Pudalov, H. Kojima, 
and M.E. Gershenson, Phys. Rev. B 75, 195330 (2007). 

53 U. Kunze, G. Lautz, Surface Sci. 142, 314 (1984). 

54 F.F. Fang, A.B. Fowler, and A. Hartstein, Phys. Rev. B 16, 4446 (1977). 

20 



55 A.B. Fowler, F. Fang, W.E. Howard, and P.J. Stiles, Phys.Rev. Lett., 16, 901 (1966). 

56 F. Fang and P.J. Stiles, Phys. Rev. 174, 823 (1968). 

57 J.L. Smith and P.J. Stiles, Phys. Rev. Lett., 29, 102 (1972). 

58 Our revision of the analysis of SdH data (described in the subsection IIII A 1|) modifies the m* 
values only at low densities (r s > 4) and does not affect m*{r s — ► 0) = 0.205m e , adopted as the 
best estimate for mj. 

59 The chosen value = 0.205m e does not contradict the cyclotron resonance data. Due to the 
Kohn theorem, the mass measured in the cyclotron resonance experiments is believed to be 
unrenormalized and equal to the band mass (and hence, independent of density). Abstreiter et 
al. [Phys. Rev. B 14, 2480 (1976)] observed the far-infrared cyclotron resonance and obtained 
m c = (0.197±0.005)m e independent of n for n > 1 x 10 12 cm~ 2 . However, for n < 1 x 10 12 cm~ 2 , 
they observed a sample-dependent noticeable increase of m c with decreasing density (up to 
0.215m e at n = 5 x 10 11 cm -2 ), which was attributed to the localization effects. 

60 J. Huang, J.S. Xia, D. C. Tsui, L.N. Pfeiffer, and K.W. West, Phys. Rev. Lett. 98, 226801 
(2007). 

61 V.M.Pudalov, G. Brunthaler, A. Prinz, and G. Bauer, Phys. Rev. B 60 R2154 (1999). 

62 For a two-valley system with A v S> T, the total number of interaction channels is 4 x 2 = 8, 
among them 1 singlet and 7 triplet terms. 

63 X.P.A. Gao, A. P. Mills, Jr., A. P. Ramirez, L.N. Pfeiffer, and K.W. West, Phys. Rev. Lett. 
88, 166803 (2002). 

64 V. M. Pudalov, M. Gershenson, and H. Kojima, "On the electron-electron interactions in two 
dimensions" , in: Fundamental Problems of Mesoscopic Physics, eds. I. V. Lerner et al. (Kluwer 
2004), pp. 309-327. 

65 P.M. Mensz and R. G. Wheeler, Phys. Rev. B 35, 2844 (1987). 

66 H. Mathur and H.U. Baranger, Phys. Rev. B 64, 235325 (2001). 

67 G.M. Minkov et al., Phys. Rev. B 70, 035304 (2004). 

68 J. S. Meyer, V. I. Fal'ko, and B. L. Altshuler, in NATO Science Series II, vol. 72, eds. I. V. Lerner 
et al. (Kluwer Academics, Dordrecht, 2002), pp. 117-164. 

69 S.V. Kravchenko et al., Solid State Comm. 116, 495-499 (2000). 



21 




< 

b 

CO 



0.6 



0.0 



4 5 



1 ' 1 1 1 — 

(b) 





1.5 



2.0 



B L -\T l ) 



2.5 



b 

CO 



0.02 - 



0.00 



It I I 


\ ( d )" 

\ 
\ 








\ 






1 


1 1 , *■» 



- 0.08 



b 

~b 

CO 



- 0.04 



0.00 



0.0 0.3 0.6 0.0 0.3 0.6 0.9 
T(K) T(K) 



FIG. 1: (color online) Shubnikov - de Haas oscillations normalized by the amplitude of the first har- 
monic A\\ (a) sample Si6-14, n = 6.1 x 10 11 cm~ 2 , T = 36 mK; (b) sample Sil-46, n = 1 x 10 12 cm~ 2 , 
T = 200 mK. Dots represent the experimental data, solid curves - the theoretical dependences (7) 
modified for a finite Ay and calculated for Ay = 0.4 K and 0.7 K for samples S16-14 and Sil- 
46 respectively. Panels (c) and (d) show the temperature dependences of the amplitude of SdH 
oscillations for Si6-14 (n = 5.5 x 10 n cm -2 ) and Sil-46 (n = 1 x 10 12 cm -2 ), solid curves - the non- 
interacting LK-model [Eqs. (JSJ)], dashed curves - the fit based on the interaction theory [Eqs. ([9l 
TUJ HI]) ]. The linearized Eq. (JT2J) is indistinguishable from Eqs. <P [TD j ITT ]) within the studied T 
range and not shown therefore. 
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FIG. 2: (color online) Temperature dependences of the conductivity cr(T, = 0) for sample Si6- 
14 at various electron densities n = 15.1,10.0,5.5,4.0,2.9,1.8, in units of 10 11 cm -2 , from top to 
bottom. Circles show experimental data, red curves - the theoretical dependences Eq. ([S} calculated 
with Ay = 0.4 K. On the left panel, the dashed blue curve corresponds to T = 0.5(1 + F$) 2 E F , 
the applicability of the ZNA theory is violated at a higher T. The right panel shows the same 
data set within a narrower temperature interval, the thin arrows correspond to T = r" 1 and 
the thick arrows - to the temperature of the crossover between ballistic and diffusive regimes, 
T* = (1 + Fq)/2itt. The dash-dotted green curves on the right panel were calculated with three 
triplet components in Eq. (1) (the valleys are completely intermixed), the blue dashed curves - 
a(T,B) = <td + Sac + -^triplet (Tt v ) x 5(Jt(T) with -/Vt r ipi ct continuously varying between 3 (for 
T < t~ 1 , Ay) and 15 (for T > r" 1 , Ay). 
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FIG. 3: (color online.) Temperature dependences of the conductivity for samples Sil-46 [(a) 
and (b)] and Si6-14 [(c)-j-(h)] in different in-plane magnetic fields (from top to bottom, Bu = 
0, 0.6, 1, 1.5, 2, 3 T in panels (a) and (b), B\\ = 0, 0.6, 1,2,3 T in panels (c)-^(h)). Experimental 
data are shown as circles, the solid curves show the theoretical dependences calculated for sample 
Si6-14 with Ay = 0.4 K and for sample Sil-46 with Ay = 0.7K. The Fq value is the only fitting 
parameter in comparison with the theory^, the corresponding values of Fq are shown in Fig. 5. 
The values of n are shown in units of 10 11 cm -2 . 
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FIG. 4: (color online.) Magnetoconductance for sample Si6-14 at different electron densities and 
temperatures. Experimental data are shown by dots, the theoretical dependences calculated ac- 
cording to Eqs. (5,6) - by solid curves. The F£ value is the only fitting parameter in comparison 
with the theory^, the corresponding values of Fq are shown in Fig. 5. Arrows indicate the fields 
corresponding to the condition g^BB\\/2Ep = 0.1. The values of n are shown in units of 10 11 cm -2 . 
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FIG. 5: (color online.) (a) The Fq values obtained from fitting the cr(T) and cr(B^) dependences 
with the theory^ (blue and red symbols, respectively). Open red squares show the shift of several 
Fq values extracted from cr(B^) if one takes into account the g(B^) dependence and suppression 
of the WL corrections by B^ (see the text). The dashed curve corresponds to Ffifrs) extracted 
from the SdH data^ using the LK theory, the dash-dotted curve - to the empirical approach 



used in Ref. l25l. The shaded regions in panels (a) and (b) show the FQ(r s ) dependence (with the 
experimental uncertainty) obtained from fitting our SdH data^ with the theory^, (b) Comparison 
of Fq values calculated from a(T, B = 0) using the same fitting procedure (see the text): - present 



work, and ▲ - data recalculated from Refs. 



36 and 



38l . respectively. 
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